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We investigate theorerically the non-equilibrium transport through a molecular quantum dot as a 
function of gate and bias voltage, taking into account the typical situation in molecular electronics. 

In this respect, our study includes asymmetries both in the capacitances and tunneling rates to the 
source and drain electrodes, as well as an infinitely large charging energy on the molecule. Our 
calculations are based on the out-of-equilibrium Non-Crossing-Approximation (NCA), which is a re¬ 
liable technique in the regime under consideration. We find that Coulomb blockade edges and Kondo 
peaks display strong renormalization in their width and intensity as a function of these asymmetries, 
and that basic expectations from Coulomb blockade theory must be taken with care in general, ex- 
pecially when Kondo physics is at play. In order to help comparison of theory to experiments, we 
also propose a simple phenomenological model which reproduces semi-quantitatively the Coulomb 
blockade edges that were numerically computed from the NCA in all regimes of parameters. 

PACS numbers: 73.63.-b, 73.63.Kv, 72.15.Qm 


I. INTRODUCTION 

In the last decade, an enormous amount of research in 
the field of nanoelectronics was devoted to the study of 
electronic transport through quantum dots. A seminal 
and very versatile type of device is the gate-defined ar¬ 
tificial atom created at the two-dimensional electron gas 
confined in semiconducting heterostructures. 1_ — More re¬ 
cently, molecular quantum dots^~— have gained an im¬ 
portant momentum, due to the possibility of playing with 
the chemistry of the molecule, which compensates the re¬ 
duced gate tunability of these systems. In all these quan¬ 
tum dot setups, several physical effects are generically 
observed when the system is cooled at cryogenic temper¬ 
atures: Coulomb blockade is probably the major one, due 
to the large charging effects in these nanoscopic quantum 
dots. In addition, more subtle correlations show up due 
to the Kondo effec t 13 ’ 14 , which implies a resonance at 
the Fermi energy in the spectral density of the dot state, 
which in turn leads to an anomalous peak in the differ¬ 
ential conductance G(V ) = dl/dV at zero bias voltage 
V. 

However, important and unexpected correlation effects 
can also take place in transport features at finite bias 
voltage. For instance, Coulomb blockade edges were 
shown by Konemann et air — to present a strong width 
renormalization in the situation of large tunnel asym¬ 
metries between source and drain electrodes, that is also 
quite typical in molecular quantum dots.— Enhancement 
of cotunneling features of excitation levels (as precursor 
of the Kondo effect) were also described by Paaske et 


aL 17 and later widely observed. All these features need to 
be incorporated in the current theoretical developments 
of non-equilibrium quantum transport, which has so far 
mainly focused in analyzing the main universal Kondo 
signatures in G(V) under many simplifying assumptions 
(such as symmetric capacitive and tunnel coupling to the 
electrodes). In particular, the experiment by Konemann 
et a/— investigated a strongly confined semiconducting 
quantum dot coupled to two conducting leads with highly 
asymmetric tunnel couplings, and studied tunneling res¬ 
onances at nearly opposite values of V. While in absence 
of interaction effects one expects nearly the same inten¬ 
sity and width of the two resonances (see Section ED, 
the Coulomb blockade edge at positive voltage was how¬ 
ever nearly five times more intense and nearly two times 
narrower than the peak at negative voltage. The au¬ 
thors assumed that for negative voltage, the width in the 
spectral density of the dot level is two times larger than 
for positive voltage. This is based on the fact that the 
charge-transfer peak (the one near the dot level E£) in 
the density of states of the Anderson model is two times 
wider when the dot is occupied . 18 ’ 19 However, as we shall 
show, the peak for negative voltage corresponds to the 
intermediate valence (IV) regime, in which the width of 
the spectral density is approximately 3/2 times larger 
than in the empty-orbital regime. In addition, the ratio 
of widths for negative and positive voltages depends also 
on the asymmetries of the source and drain capacitances. 
Our calculation will also be able to take into account the 
presence of a Kondo resonance at small bias, which in 
previous Coulomb blockade theory was clearly lacking— 
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In this paper, after a brief discussion of what is ex¬ 
pected in the absence of Coulomb interactions for asym¬ 
metric coupling to the leads, we calculate the differen¬ 
tial conductance dl/dV in the Keldysh formalism within 
the non-crossing approximation (NCA) i 20 i 21 The NCA 
reproduces well the scaling relations with temperature 
T and bias voltage in the Kondo regime^ and was 
for example also successfully used to interpret experi¬ 
mental results on a controlled crossover between SU(4) 
and SU(2) Kondo states driven by magnetic field in a 
nanoscale Si transistor—, as well as quantum phase tran¬ 
sitions involving singlet and triplet states in molecular 
quantum dots^ It was also used recently to discuss non¬ 
equilibrium Kondo spectroscopy in a system of double 
quantum dots— For asymmetric coupling to the leads, 
we will also derive an analytical expression which de¬ 
scribes semiquantitatively the NCA results. This expres¬ 
sion is used to interpret the experiments of Konemann et 
al— providing a good agreement, and could be used for 
further studies of molecular systems. 

The paper is organized as follows. In Section [TT] we 
present the Hamiltonian and the equation describing the 
current in presence of a finite bias. In Section Hill we com¬ 
pute for comparison the differential conductance G(V) in 
absence of Coulomb interaction. Section HVl contains the 
NCA results for the non-linear conductance, and presents 
an analysis of the positions and widths of the peaks 
displayed in G(V) for different experimentally relevant 
conditions. In Section El we discuss the experiment of 
Konemann et al ., and end up in Section m with a sum¬ 
mary. 


chemical potential of lead v. Note that the interchange 
of right and left leads or electrons by holes is equivalent 
to a change of sign of V. We take the arbitrary origin 
of energies at p R = 0. The capacitance effects modify 
the energy necessary to add an electron to the dot with 
the lever arm parameter a (which depends on the source, 
drain and gate capacitances) 15 i 16 , defined as: 

E d = E° d + aeU, (2) 

where 0 < a < 1, and e is the electron charge (taken 
positive). For the situation of symmetric voltage drop 
that is usually considered, a is often taken near 1/2, but 
this may not be the case in many experimental setups. 

We obtain the conductance G = dl/dV differentiating 
the current I (numerically from NCA results), which is 
given by25 

I = C J dujp(u,V,E d )[f L (uj) - 

r, 8neA L A R 

6 “ hA ’ (3) 

where /„(w) = /( to — p u ) is the Fermi distribution in each 
lead, with f(co) = l/(e u / kT + 1), and p(co,V,E d ) is the 
non-equilibrium (voltage-dependent) spectral density of 
the impurity level. 


III. NON-INTERACTING RESONANT-LEVEL 
MODEL 


II. MODEL 

We consider the Anderson impurity model with infi¬ 
nite on-site repulsion 17, an assumption that is realistic 
for most molecular quantum dots, but also well obeyed 
for very small semiconducting devices in the vicinity of 
a charge-degeneracy point. The impurity states thus in¬ 
volve a singlet configuration |s) together with a spin dou¬ 
blet \<j) (a =f or 4-) corresponding to one additional elec¬ 
tron (or hole) in the dot, taking this particle from the 
left ( L ) or right ( R ) conducting leads to which the dot is 
connected. The Hamiltonian thus reads: 

H = ^E d \a){a\ + ^e vk cl ka c vka 

<y ukcr 

+ y~!(Ifc iy | <7 )( S l C ^fccr +H.c.), (1) 

vkcr 


where the constraint |s)(s| + |cr)(cr| = 1 is imposed. 
Here cj) fc(T create conduction states in the lead v with 
spin a and wave vector k. The tunnel couplings of the 
quantum dot to the leads are A v = i r J2k |14 l/ | 2 J(u; — e v k), 
assumed independent of energy w, and we define the total 
hybridization A = Ai + A R . We take the sign of the bias 
voltage in such a way that pl~Pr = eV, where p v is the 


Before discussing in more detail the complex effects of 
Coulomb interactions, we present in this section a simpli¬ 
fied model that sheds light onto the expected behavior of 
the differential conductance G(V). The first assumption 
that we make in this section is that the tunnel couplings 
to the leads are very asymmetric. This hypothesis is 
realistic for molecular quantum dots and also in the ex¬ 
periment of Konemann et al. in Ref. [TBl . We make here 
the choice A R 3> A Rl so that the system is in equilib¬ 
rium with the right lead and the explicit dependence on 
the voltage of the spectral density p(u>) disappears [the 
dependence of p(ui) with V through E d given by Eq. © 
remains]. Then one has 

p(w, V, E d ) ~ p R (u, E d ), (4) 

where p R (uj,E d ) is the spectral density of the impurity 
states in equilibrium with the right lead, for a dot energy 
given by Eq. ©■ Using equations m, one obtains at 
T = 0 


G(V) = Cep R (eV, E d + aeV) 
Jo dE d 


( 5 ) 


In the limit of asymmetric capacitances a -A 0, only the 
first term survives and G(V) is just a map of the spectral 
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density p R (cu,E®), as is usual with experiments with a 
scanning tunneling microscope (STM). Therefore, in this 
limit, G(V) has a single peak located at V = E d /e (the 
charge-transfer peak), and no peak for the opposite bias 
voltage. In addition, a Kondo peak at V = 0 can show 
up in presence of Coulomb interactions, if —Efj A. 

The main assumption that we make in this section is 
to take a non-interacting resonant level spectral density 
of width 2A at half maximum, p R (u>, Ed) — Po(oj — Ed), 
with 


p 0 (u) - E d ) 


A/n 

(w - E d Y + A 2 ■ 


( 6 ) 


This is a crude approximation which is certainly not 
valid in presence of interactions (it misses the important 
Kondo anomaly), but it can describe qualitatively the po¬ 
sition and width of the Coulomb peaks in G(V) derived 
from the charge-transfer peak in the spectral density of 
the dot states, which can be renormalized from bare val¬ 
ues by the interaction, as we discuss in Section IIVB 21 
Assuming, as in Eq. ©, that dp R /dE d = -dp R /dui 
[but not p R = po yet], the second term of Eq. ([SJ can be 
easily evaluated and one obtains 


G{V) ~ Ce[(l -a)p R (eV,E d ) +ap R (0,E d )}, (7) 



FIG. 1: (Color online) Equilibrium (V = 0) spectral density of 
the localized state for the infinite-interaction Anderson model 
in the empty orbital regime (dashed line: E d = 4A, T = 
0.01 A), in the intermediate valence (IV) regime (dashed-dot 
line: E d = —0.733A, T = A),— and in the local moment 
regime (full line, E d = — 4A, T = 0.01A) that shows an 
additional Kondo feature at low energy. 


noting that E d displays an explicit dependance upon the 
voltage from Eq. ©. While the assumption made above 
for the derivatives of the current is not exactly true (even 
in equilibrium) because of the dependence of the weight 
and width of the peaks with E d , this simple expression 
reflects a general feature of the expected Coulomb peaks, 
in connection to the charge-transfer peak of the equilib¬ 
rium spectral density p R . Indeed, now G(V) displays two 
finite bias peaks, one that is similar to the STM case and 
is located at voltage eV = E d [namely eV = E d /(l — a)\, 
and a second “intermediate valence” (IV) peak at an op¬ 
posite voltage value such that E d = E® + aeV = 0. In 
this latter case, the impurity experiences indeed maximal 
charge fluctuations because the level is located exactly 
on-resonance, and thus one can expect a different width 
renormalization than for the STM peak, due to the dif¬ 
ferent influence of the Coulomb interaction, as we will 
demonstrate later on. For positive E d (corresponding 
to a vanishing charge in the dot, also called the empty 
orbital regime), the STM peak lies at positive bias and 
the IV peak is located at negative bias, and conversely 
for negative E° (corresponding to a spin-degenerate odd 
charge in the dot, also called the local moment regime). 
Thus, if p R = po is assumed (which is valid only in ab¬ 
sence of Coulomb interaction), the conductance G(V) 
displays two peaks as a function of voltage: i) an “STM- 
like” peak at V = E®/[( 1 — a)e) of width 2A/[(1 — a)e] 
and maximum intensity (1 — a)eC/ A; ii) an IV peak 
at V = —E d /(<xe) of width 2A/(ae) and maximum in¬ 
tensity aeC/A. Therefore, in the absence of Coulomb 
interactions, the amplitude and width renormalizations 
are related, and cannot reproduce the experimental ob¬ 


servations made in Ref. fl5l . We now turn to a full mi¬ 
croscopic calculation of the problem, without recourse to 
the many assumptions made above, that incorporates the 
correct nature of the Coulomb blockade edges and Kondo 
resonances for all parameters in the device. 

IV. NCA RESULTS FOR THE INTERACTING 
MODEL 

A. The equilibrium spectral density 

We consider now Hamiltonian ([1]) in its full complex¬ 
ity for the case of a finite bias voltage. The non-crossing 
approximation (NCA) has proved in the past an excel¬ 
lent technique to deal with large Coulomb interactions, 
and can be extended to the present non-equilibrium sit¬ 
uation. In order to gain intuition, we first focus on the 
zero-bias (V = 0) spectral density of the localized level, 
which is shown on Fig. [T] for low temperatures and three 
values of the dot energy level E d , corresponding to the 
empty orbital ( E d = 4A), local moment ( E d = —4A), 
and mixed valence ( E d ~ 0) regimes^ For positive E d 
the result is very similar to the non-interacting resonant 
level prediction, Eq. ©, with a spectral width 2A, ex¬ 
cept for the fact that the position of the peak is slightly 
shifted to higher energy Ey = 4.72A. This energy shift 
E d B — E d is in very good agreement with the Haldane 
shifts^! (A/7r)ln(D/A) = 0.733A for the half-band width 
D = 10A that we have used in the calculations. This 
is expected because the NCA incorporates naturally the 














lowest order perturbative terms of the development in 
the hybridization V%. Note that the small spike at the 
Fermi energy w = 0 is an artifact of the NCA in the 
empty orbital regime. 

In the local moment regime at negative E the charge- 
transfer peak is nearly two times wider due to the un¬ 
derlying extra spin degeneracy, as was discussed in the 
introduction and in previous work s 18 ’ 19 , and the NCA 
correctly reproduces this behavior. The total intensity 
is also reduced by correlation effects, as discussed phe¬ 
nomenologically in Section IIV B 21 and the position of 
the peak is also shifted to higher energies. More impor¬ 
tantly, the screening of the spin-degenerate states leads 
to a Kondo peak at the Fermi level, a non-trivial effect 
that the NCA is able to describe properly. 

Finally, the density of states for the mixed valence sit¬ 
uation (E^ s = 0) shows a peak centered at the Fermi 
level as expected, with a width of about 3A that inter¬ 
polates between the widths found in the empty orbital 
and local moment regimes. The temperature in this cal¬ 
culation has been increased to avoid a shortcoming of the 
NCA^ 

Regarding the conductance G(V), these results antic¬ 
ipate similarities and differences with the naive expecta¬ 
tions based on the non-interacting resonant level model 
discussed at the end of the previous section (apart from 
the obvious lack of Kondo resonance in the latter model). 
A minor point is that E® should be replaced by an effec¬ 
tive level E®’ e , which includes the Haldane shift. Based 
on Fig. [I] and Eq. (JTJ) the width of the STM-like peak in 
the empty orbital regime should be 2A/[(1—a)e], thus 4A 
at a = 1/2 for a symmetric voltage drop (as for the non¬ 
interacting resonant level model). However, the width of 
the STM peak in the local moment regime should rather 
be 4A/[(1 — a)e], thus 8A at a = 1/2. More subtly, the 
IV peak (located at opposite voltage to the STM peak) 
should have the same width for the empty orbital and 
local moment regimes, namely 3A/[(1 — a)e], thus 6A 
at a = 1/2. These general trends, that we will now ex¬ 
amine microscopically from full NCA calculations of the 
conductance, show that the experimental analysis of the 
width of Coulomb blockade peak must be done with care. 


B. Conductance for asymmetric tunnel coupling to 
the leads 

1. Numerical NCA results 

In this subsection, we calculate the conductance using 
the non-equilibrium formalism within the non-crossing 
approximation (NCA) 20 ’ 21 for a large choice of parame¬ 
ters. Details on the technique are given for example in 
Ref. jH- We choose here asymmetric tunnel coupling to 
the leads Ar = 60 Ar, and set the reference energy scale 
as A = Ar + Ar = 1, with a temperature T = 0.01 and a 
half-bandwidth of the leads D = 10. Owing to the large 
ratio Ar/Ar 1, the impurity is effectively in equilib¬ 



FIG. 2: (Color online) Finite bias differential conductance 
G(V) in the empty orbital regime E° = 4A = 4 for several val¬ 
ues of the capacitance ratio a = 1/2,1/3,1/4 defined accord- 
in □ 

(t 



FIG. 3: (Color online) Finite bias differential conductance 
G(V) in the local moment regime E° = —4A = —4 for 
several values of the capacitance ratio a = 1/2,1/3,1/4 de¬ 
fined according to Eq. ©, and with an hybridization ratio 
Ar/Al = 60. 


rium with the right lead, although the calculations were 
in practice performed using the full out-of-equilibrium 
NCA. 

In Fig. [2] we show G(V) in units of the conductance 
quantum Go = 2 e 2 /h for positive (empty orbital 
regime) and several values of capacitance asymmetries 
a. The small spike at V = 0 is again an artifact of 
the NCA and should be disregarded. For symmetric 
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voltage drop, a = 1/2, the discussion of Section Hill 
based on the non-interacting resonant level predicts two 
peaks of equal width 4A and equal height eC/ (2A) at 
eV = ±2 E° d ' cS « 9.4. However, the more refined analysis 
of Section IIV A1 based on the general features of the in¬ 
teracting spectral density, predicts different width renor¬ 
malization: a width 4A for the STM peak, and a width 
6A for the IV peak. These trends are well described by 
the full NCA calculation shown in Fig. [2] However, the 
conductance G = dl/dV is not directly given by the spec¬ 
tral density, but the dependences of its effective width 
and weight with E d affect its shape [second term of Eq. 
(5)]. This is particularly important for the IV peak, for 
which these dependences are larger. This will be dis- 
ussed in more detail in the next subsection. The width 
of the IV peak within the NCA is about 20% larger than 
6A. This can be due to the above mentioned effects or 
to the limitations of the NCA at low temperatures in 
the IV regime. In any case, our results demonstrate that 
the STM and IV peaks are generically asymmetrical, and 
that their widths depend sensitively on the charge state 
present on the quantum dot. 

As a decreases, the tendencies predicted by the simple 
analysis of Section Hill remain qualitatively valid. The 
position of the STM peak tends to and its width to 
2A in the limit a —> 0. In contrast, the IV peak broadens 
enormously and moves to more negative values of V when 
the capacitances become very asymmetric. 

In Fig. [3] we show the results for the local moment 
regime at negative E®. In this case, both peaks are much 
broader than the predictions of the naive expectations 
based on a non-interacting resonant-level spectral den¬ 
sity. Indeed, the STM peak now lying at negative volt¬ 
age shows strong interaction effects (already witnessed in 
the spectral density in Fig. [1] that make it twice broader 
than without interactions (the width is 8A for symmet¬ 
ric voltage drop a = 1/2). The IV peak, now located at 
positive bias voltage, displays again a width that is about 
3/2 wider than the simple expectations outlined at the 
end of Section m due to combined charge fluctuations 
and Coulomb effects (the width is 6A for a = 1/2). The 
IV peak is also shifted to the left, since for a = 1/2 one 
would expect it to lie at eV = 2 E®’ eS « 6.5A, while it is 
located at 5.7A. While naively one would expect a ratio 
of maxima inversely proportional to the width, the inten¬ 
sity of the STM peak is smaller than expected, especially 
for a = 1/2. This further correlation effect is discussed 
in the next subsection where we present a more refined 
analytical approach. Again, as a decreases, the width 
and intensity of both peaks evolve in agreement with the 
expectations of the non-interacting resonant-level model. 


2. Phenomenological approach 


The main reason for the failure of the non-interacting 
resonant-level model to account for Coulomb blockade 
peaks at asymmetric tunnel coupling to the leads is that 
the correct width and weight of the peak in the impurity 
level spectral density p(ui) both depend in a non-trivial 
way on the occupancy, as discussed in Section HV Al One 
of the simplest ways to modelize this effect is the alloy 
analogy approach to periodic Anderson models, in which 
the electrons with spin down are frozen to calculate the 
dynamics of electrons with spin up and vice versa i 29 ’ 30 
This approach is also related to the so-called Hubbard- 
III approximation*^! in which the weight of the charge- 
transfer peak near Ed for spin a (or the lower Hubbard 
band in a Hubbard model) is proportional to the number 
of unoccupied states with the opposite spin, 1 — (n_ CT ). 

To correct the width inconsistency in the non¬ 
interacting resonant level model, we have fitted the width 
of the NCA peak as a function of Ed at a temperature 
T = 2A (such that the Kondo peak is absent and the 
width of the IV peak is well behaved^), with other pa¬ 
rameters taken as in Section IIV B II For simplicity, in 
this section and in section El we replace the notation 
E^ s by Ed , incorporating the Haldane shift in the lat¬ 
ter. From this, we obtain that the effective half-width at 
half-maximum can be approximated as: 


A e ff 

ET 


3 

2 


-tanh 

2 


E d \ 

2.77A ) 


( 8 ) 


This simple formula nicely interpolates between the 
empty orbital case (A e ff = A for Ed A) and the lo¬ 
cal moment regime (A e ff = 2A for —Ed —A), passing 
through the mixed valence regime (A e g = (3/2)A for 
\E d \ < A). 

The above discussion suggests that a better description 
of the single peak in the equilibrium local density of state 
can be obtained using the improved formula: 


Pr{w, E d ) 


(1 - (n a ))A eS /TT 
(ui — Ed) 2 + A“ ff 


(9) 


where (n a ) = f dujpji(uj)fR(u}) is the occupation of the 
dot for a given spin, independent of the spin in the ab¬ 
sence of an applied magnetic field. Integration gives (n a ) 
in terms of the digamma function ip(x): 


(n>a) 

4V 


1/2-ttfl 

3/2- 

—Imip 

7r 


A e ff T i(E d Pis) 
2 ttT 


( 10 ) 


As expected, the most drastic difference with the non¬ 
interacting resonant-level model for negative E® is the 
appearance of the Kondo peak at V = 0. Note that 
its main features are independent of a, which will be 
discussed in Section IIV C 21 We now turn to a physical 
interpretation of the numerical NCA result. 


Replacing then Eq. J9| in Eq. ([3]) gives the current: 

I = C( 1 -<«,»(**-*£), (11) 

Finally differentiating this expression, one obtains an an¬ 
alytical although lengthy expression for the differential 


















FIG. 4: (Color online) Differential conductance in the empty 
orbital regime E° = 4A = 4, for asymmetric tunnel couplings 
(An = 60 A_l), low temperature T = 0.01A, and two val¬ 
ues of the capacitance ratio a = 1/2,1/4. Full lines are the 
NCA simulations, and dashed lines are the results from the 
phenomenological approach. 


conductance G(V), which we do not reproduce here. At 
T = 0, G(V) contains two terms which have the form 
of those included in Eq. 0, plus some additional terms 
due to the dependence of A c ff and ( rida ) with respect to 
Ed- These additional terms affect the IV peak, for which 
the variation of A c ff and (rida) is significant. 

In Fig. [4] we compare the conductance G(V) obtained 
through this phenomenological approach to the full NCA 
result in the empty orbital case E® = 4A (note that the 
Haldane shift is included by hand in the level position 
within the phenomenological approach). We remind the 
reader that for a = 1/2, the non-interacting resonant- 
level model predicts a symmetric G(V) = G(—V) con¬ 
ductance curve, which turns out to be incorrect in 
the presence of strong Coulomb interaction. The phe¬ 
nomenological approach is thus a considerable improve¬ 
ment for the IV peak at negative voltage, although some 
quantitative discrepancies remain. In particular, the ana¬ 
lytical expression gives an intermediate valence peak a lit¬ 
tle bit narrower and shifted towards V = 0, but the gen¬ 
eral trends are well reproduced by the simple approach. 
In addition, this discrepancy might be due to limitations 
of the NCA at low temperatures in this regime. 

In the local moment regime for negative E as shown 
in Fig. 0 the comparison is less satisfactory. The Kondo 
peak is obviously missing in the phenomenological ap¬ 
proach and the positions of the Coulomb peaks show 
also some deviations, particularly for the IV peak now 
located at positive V. Nevertheless, the intensities and 
the widths (which are larger than those predicted by the 
non-interacting resonant-level model) are semiquantita- 
tively reproduced. The small discrepancy in the position 


FIG. 5: (Color online) Differential conductance in the local 
moment regime E ° = —4A = —4, with the same parameters 
as in Fig. [4] Full lines are the NCA simulations, and dashed 
lines are the results from the phenomenological approach. 


of the STM peak at negative V might be due to some 
dependence with Ed of the Haldane shift, being smaller 
for negative Ed- Having clarified the physics at play in 
the asymmetric situation, which simplifies because the 
system is in equilibrium with a single lead, we now con¬ 
sider the problem in is full complexity for the regime of 
symmetric tunnel couplings. 


C. Conductance for symmetric tunnel coupling to 
the leads 

1. Numerical NCA results 

For symmetric tunnel coupling to the leads, an ap¬ 
proximation like Eq. 0, which would permit to reduce 
the problem to the calculation of the spectral density at 
equilibrium, is no more possible. In Fig. [6] we show the 
results of the non-equilibrium NCA calculation for sym¬ 
metric tunnel barriers, A l = A r = A/2, and parameters 
similar to those of the experiment by Konemann et al 
but with negative E In contrast to the results of Section 
IIV B 11 the symmetry now imposes that G(V) = G(—V) 
for Al = Ar and a = 1/2. 

While one would expect that the non-interacting 
resonant-level model is meaningless in this non¬ 
equilibrium interacting case, the symmetric coupling to 
the leads “blurs” in the conductance the asymmetry due 
to correlation effects. The simplest non-interacting reso¬ 
nant level model with renormalized parameters is indeed 
able to reproduce the Coulomb blockade peaks obtained 
with the NCA (of course the Kondo anomaly is not cap¬ 
tured). Except for the small modification of the a value, 
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FIG. 6: (Color online) Differential conductance at symmetric 
tunnel couplings, A l = A_r, for parameters close to the exper¬ 
imental measurement of Konemann et ah —, namely A = 1.15 
meV, E° = —4.5 meV, T = 0.014 meV and a = 0.425. 
The full line is the non-equilibrium NCA calculation, while 
the dashed line corresponds to a fit to the non-interacting 
resonant-level model with renormalized parameters A cfi =2.1 
meV, E® = —3.44 meV, a — 0.45 and a scale reduction factor 
0.355. 

which is hard to justify, the changes in the other param¬ 
eters of the model are expected: a width 4.2 meV nearly 
two times larger than expected (2.30 meV) because of 
correlation effects, an effective Ed shifted upwards due to 
renormalization effects (the Haldane shifts is 1.04 meV 
for the half band width D = 20 meV used in the fig¬ 
ure) and smaller intensity also due to correlation effects 
(weight proportional to 1 — (n- c ,)). We now turn to a 
detailed description of the Kondo resonance for a large 
range of parameters. 


2. On the width of the non-equilibrium Kondo peak 

In this section, we study the relation between the width 
of the low bias Kondo resonance in the differential con¬ 
ductance G(V) (which is present for negative enough E°) 
and the Kondo peak in the equilibrium spectral density 
of the impurity level p(u>). We define Tr (resp. T£) 
as the energy scale corresponding to half the width at 
half maximum of the peak in G(V) (resp. p(co)) near 
V = 0 (resp. u> = 0). Both quantities are of the or¬ 
der of the Kondo temperature, but differ from a non¬ 
uni versal prefactor. Since the Kondo temperature is of¬ 
ten extracted experimentally from the differential con¬ 
ductance curve, but theoretically the equilibrium density 
of states is more standardly computed, for instance with 
the Numerical Renormalization Group (which has diffi¬ 
culties in grasping the effect of finite bias voltage), we 


TABLE I: Ratio Tr/Tr of the Kondo temperatures extracted 
from finite bias conductance G(V) and equilibrium density of 
states p(cu), obtained for different values of the level energy 
Ed and various asymmetry ratios Ah/A l- 


Ar/Al 

E d 

-2 

-3 

-4 

-5 

-6 

1 

1.59 

1.45 

1.33 

1.31 

1.30 

2 

1.50 

1.19 

1.078 

1.09 

1.08 

5 

1.29 

1.03 

1.008 

1.003 

- 1 

10 

1.29 

~ 1 

~ 1 

~ 1 

- 1 



FIG. 7: (Color online) Differential conductance near G(V) 
(broken lines in red, corresponding to the left vertical scale), 
and equilibrim spectral density of the quantum dot p{oi) (full 
line in black, corresponding to the right vertical scale), com¬ 
puted for A = 1, Ed = —4 and several tunnel asymmetry ra¬ 
tios R = Ah/A l = 1,5,30. The coefficient A = 4AhA_l/A 2 
is introduced to normalize the conductance to the density of 
states. 

consider here systematically the ratio Tr/Tr, for sev¬ 
eral ratios of couplings to the leads Ar/Al, keeping 
A = Ar + Al = 1 fixed, and also for several values of the 
level position Ed- In this section we keep Ed fixed and 
assume a symmetric voltage drop pl = ~Pr = eV/2, as 
V is increased, and take the half band width D = 10. 
Both the differential conductance G(V) and equilibrium 
density of states p(co) are computed with the NCA, the 
former with the non-equilibrium formalism, and the lat¬ 
ter at equilibriumi 20 ! 21 ’ 28 The resulting ratio of Kondo 
temperatures is shown in Table (U 

It can be noticed from Table [T] that the ratio Tr/Tr > 
1 and decreases towards unity when both the tunnel 

















asymmetry ratio R = Ar/Al is increased and the im¬ 
purity level position E d is decreased negatively, which 
can be understood as follows. For large asymmetry of 
the tunnel coupling to the leads, R^> 1, the dot can be 
considered at equilibrium with the right lead, and Eq. 
0 becomes valid. Furthermore, in the Kondo limit for 
fixed chemical potential p, the Kondo peak is insensitive 
to a small change in the localized level E d , and remains 
close to the Fermi level. For |eV| <C this is supported 
by Fermi-liquid approaches ! 32 ’ 33 Thus, one can drop the 
second term in Eq. ([5]) and G(V) for |eVj < gives 
directly the spectral density of the dot state pR^eV, E d ). 
This effect is evident in Fig. [3 where the evolution of the 
Kondo peak in conductance G(V) with increasing asym¬ 
metry ratio R is shown and compared to the equilibrium 
density of states pR(eV,E d ). 


V. COMPARISON WITH EXPERIMENT 

In this last section, we discuss the experiment of 
Konemann et This experiment corresponds to a 

value of E d that is positive and very large (two orders of 
magnitude bigger than A). This situation is very hard for 
the numerical solution of the self-consistent NCA equa¬ 
tions, but since the Kondo effect does not take place, 
we can compare the experiment to the phenomenolog¬ 
ical approach developed in Section IIVB21 which takes 
into account the changes in both width and height of the 
charge-transfer peak in the spectral density with the dot 
level E d - The experiment was done for highly asymmetric 
tunnel coupling to the leads, so that the main assumption 
of Section IIVB 21 is fulfilled. 

In the experiment, there is a peak at V =-9.7 mV with 
a half width of 76 pV and another one at V =8.9 mV 
with a half width of 35.5 pV which is about five times 
more intense. For a cs 1/2 (as suggested in Ref. [l5h . 
the non-interacting resonant level model incorrectly pre¬ 
dicts two peaks of equal width and height at opposite 
voltages (see the end of Section m- For a = 1/2, our 
more accurate phenomelogical theory predicts a ratio of 
about 3/2 between the width of the IV peak at nega¬ 
tive V and the STM peak at positive V , and not a ratio 
of about two. However, the fact that the positions of 
the two peaks are not symmetric with respect to zero 
bias shows that a is in fact different from 1/2. Using 
a = 0.45 to adjust the correct peak position, this pro¬ 
vides an additional contribution to the ratio in the peak 
widths. In Fig. [5] we display our results for parameters 
close to the experiment. While we are unable to repro¬ 
duce exactly the intensities, positions and widths of both 
peaks, we obtain a semiquantitative agreement with ex¬ 
periment. Our conductance curves can be fitted with the 
sum of two Lorentzians: one located at V = —9.97 mV 
with half-width at half-maximum 64 pV and peak value 
1.17 A, another one located at V = 8.18 mV, with half¬ 
width at half-maximum 36 pV and peak value 4.21 A. In 
the future, comparisons to further experimental data for 



V(meV) 


FIG. 8: (Color online) Differential conductance using the phe¬ 
nomenological model Eq. CD for parameters corresponding 
to the experiment by Konemann et al.<^ A = 0.02 meV, 
E d = 4.5 meV, T — 0.01 meV and a = 0.45. Full (dotted) 
line corresponds to negative (positive) bias voltage. 


molecular quantum dots deep in the Kondo regime would 
be interesting, in order to test our predictions. 


VI. SUMMARY 

Motivated by experiments with semiconducting and 
molecular quantum dot o 15 ’ 16 , we have analyzed the com¬ 
plete shape of the differential conductance G{V) in the 
infinite-!/ Anderson model connected to two leads, as a 
function of the ratio of the potential drops between the 
dot and both leads, and for various tunneling asymme¬ 
tries. The simplest situation is the one found in scanning 
tunneling spectroscopies, where both the potential drop 
and tunneling rates are strongly asymmetric towards the 
metallic surface on which a molecule is deposited. As is 
well-known in this situation, the conductance provides 
exactly a measure of the equilibrium spectral density of 
the quantum dot and displays only one Coulomb peak 
when voltage is aligned to the charge-transfer peak in the 
spectral density (this peak position may contain renor¬ 
malization effects). When the impurity level is deeply 
negative, a local moment forms, and the Kondo effect 
leads to an additional anomaly near zero bias, both in 
the density of states and differential conductance. 

The situation of an arbitrary ratio for the capacitance 
to the source and drain is more complex (even still as¬ 
suming very asymmetric tunneling rates). We find that 
the Coulomb peak position and width in this regime not 
only display trivial changes expected from the modifi- 
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cation of the capacitance ratio, but are also corrected 
by many-body effects. As it is known from previous 
calculations , 18 i 19 , and reproduced by our NCA results, 
the width of the charge-transfer peak is about two times 
larger for large and negative level position (local moment 
regime) than for large and positive position (empty or¬ 
bital regime). In addition, an extra Coulomb peak ap¬ 
pears at opposite voltage values to the peak found in 
the STM situation, as expected from Coulomb blockade 
theory. However, this peak has in general a different na¬ 
ture, and relates to the intermediate valence situation of 
the underlying quantum impurity model, and displays ac¬ 
cordingly a width renormalization of about 3/2 compared 
to expectations based on a non-interacting effective res¬ 
onant level model. All these effects can be taken into ac¬ 
count by a simple phenomenological approach described 
in Section IIVB 21 which is able to describe well our nu¬ 
merical simulations (except for the presence of the Kondo 
resonance) and agrees semiquantitatively with some pub¬ 
lished experimental results [see Section IV]. 

We also provided some general study of the Kondo 


temperature extracted from the differential conductance 
(as typically done in experiments), and made some pre¬ 
cise connection to the Kondo scale extracted theoretically 
from the equilibrium spectral density. We find that the 
former scale can be up to 60% larger than the latter one, 
and that they coincide only for asymmetric tunnel cou¬ 
plings deeply in the Kondo regime. 

We hope that all these results will serve as guidelines 
for interpreting more quantitatively experimental data 
obtained from molecular quantum dots, since many-body 
effects can quantitatively alter the predictions of basic 
Coulomb blockade theory. 
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